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Abstract. We present the constraints on inflationary parameters in a flat 
ACDM universe obtained by WMAP three year data release, plus smaller scale 
CMB and two LSS data sets, 2dF and SDSS (treated separately). We use 
a Markov Chain Monte Carlo (MCMC) technique combined with an analytic 
description of the inflationary spectra in terms of the horizon flow functions 
(HFF). By imposing a consistency condition for the tensor-to-scalar ratio, we 
study the constraints both on single field standard inflation and on inflation with 
the violation of the null energy condition, which leads to a blue spectrum for 
gravitational waves. For standard inflation, the constraint on the tensor-to-scalar 
ratio we obtain from CMB data and 2dF05 is: ro.oi < 0.26 at 2 cr cl. Without 
the consistency condition between the tensor-to-scalar ratio and the tensor slope, 
the constraints on the tensor amplitude is not significantly changed, but the 
constraints on the HFFs are significantly relaxed. We then show that when the 
third HFF ea is allowed to be non-zero and to be of order unity, a large negative 
(at 2(t) value for the running of the scalar spectral index in standard inflation is 
found in any set of data we consider. 



PACS numbers: CMBR theory, inflation 
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Introduction 

Inflation is the most promising paradigm for understanding the initial conditions 
for cosmic structure formation and the pattern of anisotropics in temperature and 
polarization of the cosmic microwave background (CMB). Its predictions for the 
simplest case of a single scalar field model with a nearly scale-invariant spectrum 
of curvature and tensor perturbations are in good agreement with most of the 
observational data. 

Among observations, data on CMB anisotropics have been the most selective for 
inflation, starting from COBE-DMR [T] to the NASA mission presently in operation 
WMAP . The WMAP first-year data [21 S] complemented with smaller scale data 
started to discriminate among the inflationary models and to constrain the possibility 
that the spectra may have not a pure power-law shape [SJ [6] . 

In this paper we compare inflationary predictions with WMAP three year data 
release O HI [9] plus smaller scale CMB data (which we shall denote collectively as 
CMBsmall), such as CBI [TO], ACBAR [H], VSA [H] and B2K [HUKTl]. We shall 
also use separately two LSS data sets, SDSS and 2dF [iTlfTS] . 

We compare inflationary predictions with observations, by adopting a 
parametrization of the primordial power spectra (PS henceforth) of curvature and 
tensor perturbations as: 



where ? = S , T stands for scalar and tensor, respectively, /c* is the pivot scale. We shall 
use Ai^rii, ai in terms of of the Hubble parameter H and the horizon flow functions Ci 
(HFF henceforth) evaluated at the pivot scale fc*. At the present state of the art, the 
PS of primordial perturbations in most of the single field inflationary models depend 
only on the physics at the Hubble crossing since adiabaticity on large scale is preserved. 

The HFFs are defined as ei = -H/H^ and e^+i = e,/{Hei) = [du/dN)/ti with 
i>l and N the number fo e- folds {dN — Hdt) j^Qj. For a Klein-Gordon scalar field 
(j), H and ei are related to the potential and its derivatives as: 



where Mpi = (STrG)"^/^ ~ 2.4 x 10^^ Gev is the reduced Planck mass. By comparing 
constraints on H ,ei derived by observations, the shape of the inflationary potential 
can be constrained through Eqs. ([2]). 

We use the most advanced analytic descriptions of inflationary spectra which 
relate {Ai,ni, ai) to {H,ej). These expressions are obtained by the Green's function 
method (henceforth GFM) [3T1 [22] and by the method of comparison equations 
(henceforth MCE) [23]. These two methods provide a description of cosmological 
perturbations spectra up to second order in the HFF, allowing a good accuracy also 
in inflationary models where the HFF are not so small or constant in time, i.e. the 
slow-roll condition is not well satisfied [22l [23l |24] . Note that the GFM and MCE 
methods predict the same PSs at 0(ej), but differences of percent order arise in the 
coefficients of 0(e|) terms [23] . 




(1) 




(2) 
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In this paper we address several points of current interest in comparing 
inflationary preditions with observations. At first order in HFFs we update the 
analysis of Leach and Liddle [H to WMAP3 [3 [SI [3] and to 2dF05 [IS]. At second order 
in HFFs, we first show how the large running found by the WMAP team j^Hl [S] [S] 
is recovered by using the PS parametrization in terms of the HFFs. We also give a 
model-independent constraint on the tensor amplitude by considering the tensor-to- 
scalar ratio as a free parameter and not given in terms of the tensor spectral index. 
We then study the constraints on inflationary models which violate the null energy 
condition (NEC) and predict a blue spectrum for gravitational waves |2S]. 

Our paper is organized as follows. In section II we present our methodology 
and in section III the constraints on cosmological parameters by restricting the 
parametrization for the spectra at first order in HFFs for standard infiationary models. 
In section HI we relax the assumption that the inflaton is a Klein-Gordon (KG 
henceforth) scalar fleld and we leave the tensor-to-scalar ratio as a free parameter. 
In section IV we compare the constraints on inflationary models with a blue spectrum 
for gravitational waves with standard inflationary models. In section V we present the 
results obtained allowing running(s) of the spectral indices in standard inflationary 
models, i.e. taking into account spectra to second order in HFFs. In section VI we 
conclude. 

1. Methodology 

We use the publicly available Markov Chain Monte Carlo (MCMC) package [53] 
which makes use of the Einstein-Botzmann code CAMB [57] for the computation of 
theoretical CMB anisotropics PSs and transfer functions. We restrict ourselves to flat 
spatial sections and to three species of massless neutrinos. We use as initial conditions 
for cosmological fluctuations the growing adiabatic mode. We include the CMBsmall 
data set as by default implemented in COSMOMC. Note that small scale data by the 
same experiments have been used by the WMAP team in Ref. [9] , although there may 
be additional cross-correlations between the data sets which we consider and WMAP3 
data. We restrict our analysis to linear transfer functions and therefore the LSS data 
are automatically truncated by the MCMC to kcut/h ~ 0.15 (0.2) for 2dF (SDSS). 
The convergence diagnostic on runs with multiple chains is the R statistic by Gelman 
and Rubin implemented in the MCMC. 

We use the ln(/c//c*) expansion of the logarithm of the parametrizations in Eqs. 
([!]), in order to avoid that the spectrum becomes negative for some values of the 
inflationary parameters explored by MCMC [22]. According to Eqs. ([T]), the ln(A;/A;*) 
expansion of the power spectra is: 




(3) 



and we use both the GFM ,21] and MCE [22] to relate the coeflicients 5s to the HFFs. 
The coefficients for the scalar spectrum are: 




(4) 
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bsi=ns-l^ -2ei - e2 - 2el - {2 X + 3) - X ^ (5) 

6s2 = as = -2eie2 - (6) 
and for tensors are: 

6to- -2{C+l)ei + (-2C+4-7) 4 



-X^-2X + f2 -2 + ATojeie2, (7) 

6ti =«T - -2ei -2e2 -2 (X + 1) 6162, (8) 
6t2 = ax ^ -2eie2, (9) 

where C = In2 + 7E — 2 « —0.7296 (7e is the Euler-Mascheroni constant) and 
de2/dN = 6263. At second order the next-to-leading coefficients depend on the 
method of approximation used to study cosmological perturbations during inflation: 
X = C and Ago = Ato = are obtained within the GFM ^Tl 1^: X = D = 
i - hr3 w -0.7652, Ago = {D - C)iD + In 2) - 1/18, Ato = 2I?(-D - C) - 1/9, 
within the MCE 23J. We therefore implement the above expansion in a modified 
version of the inflationary routines developed by Samuel Leach and public at 
[http://astr onomy.sussex. ac.uk/~sleach/inflation/cambinflation. html 

As pivot scale we choose k^, = 0.6T Mpc~^ as in [6J (see however |28]); this choice 
differs from those by the WMAP team (A:* = 0.002 Mpc"^ [5ii9j and /c* = 0.05 Mpc"^ 
and is better placed at the center of the range of wavenumbers probed by 
cosmological perturbations studied in the linear regime. Although the CMB spectrum 
includes at least a convolution of the power-spectrm with spherical Bessel functions, 
an effective pivot corresponding to fc» can be defined [301 121] ■ 

^ 3 ^ ^ ^ -1 (^0) 

2 1 + 0.0841nf7„i * V ; 

Since £^ ~ 145 corresponding to fc, — 0.01 Mpc~^ does not represent an interesting 
reference scale to measure the tensor-to-scalar ratio, we shall also quote in our results 
both rk, and the ratio of tensor to scalar contribution to Ce and Z = 10, denoted by 
-Rio (= CIq/Ciq). 

As parameter describing reionization we shall vary the redshift z^e, which is related 
to the optical depth r by an integration over the redshift z |26j 

nejz) 
/o "■'Hiz){l + zy 

with aT as the Thompson cross section and n,, is the electron density number. As for 
the tensor-to-scalar ratios r and i?io, we also give the constraints on r as a derived 
parameter. 

2. First Order HFF Results 

In this section we explore the constraints for a model in which the infiaton is described 
by a Klein-Gordon scalar field and the spectra in Eqs. ([3]) are taken to first order in 
HFFs (i.e. 632 = 6t2 = 0). When the infiaton is described by a Klein-Gordon scalar 
field the tensor-to-scalar ratio is not a free parameter, but is given by the so-called 
consistency condition which is: 

^ = 16ei. (12) 



I rf^TTTTT^f^, (11) 
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In this case we vary 7 free parameters: 

{Cih h\ ncdmh^Ho, 10^° As, Zre, 61, 63} , (13) 

where fib , f^crfm are baryon and cold dark matter densities, Hq is the present Hubble 
parameter, h ~ iJo/(100kms~^ Mpc^^), As(fc,) is the amplitude of super-Hubble 
curvature perturbation deep in the radiation era, ei , £2 are the first two HFFs. 
The prior on these parameters are: 0.005 < V.^h'^ < 0.1, 0.01 < ricdm < 0.99, 
40 < iJo/kms"^Mpc"^ < 100, A < Zre < 30 as in [.6j , lO"'^ < ci < 0.1; 
-0.25 < £2 < 0.2 

The one dimensional posterior marginalized probability distributions obtained 
from a Markov chain of (at least) 10^ elements (which uses the covariance matrix of 
a preliminary chain) for sampled and derived parameters are shown in Fig. 1,2 and 
in Table 1. We compare the results obtained by using the first year and the three 
years of WMAP data (plus CMBsmaU and 2dF02 [T7'), as well as the results of the 
three years of WMAP data and CMBsmaU with different LSS data sets (2dF02 [17], 
SDSS [E] and 2dF05 [E]). Different estimates of cosmological parameters obtained 
with different LSS data sets are very known (see Table 6 of [9]) and may be due to 
different selection of the SDSS and 2dF surveys: this is the reason why we do not 
combine different LSS data sets. 

One of the main results of WMAP3, i.e. a better estimate (and different from 
WMAPl) of reionization, has changed the value of the scalar spectral index ns- 
By combining with CMBsmaU plus 2dF02, WMAP3 fully marginalized 2a value is 
ns = 0.960±;]:[J^f (compared to the WMAPl correspondent value ns = 0.979±[]:[J|^). 
As is clear from the right panel of Figs. 3,4, without marginalizing on the tensor- 
to-scalar ratio, a pure Harrison-Zeldovich spectrum for scalar perturbations with no 
tensors {ug = 1 , r = 0) is disfavoured at more than 2cr by the combination of CMB 
data with 2dF02 or SDSS, but is within 2a for CMB data plus 2dF05. Another 
interesting result we find is that inflationary models predicting ns = 1 (£2 = ~2£i) 
and tensors are also disfavoured at 2a by the same data sets for which the HZ spectrum 
is disfavoured with the same cl (again, this is not true for CMB data plus 2dF05). 
As shown by Fig. 4, this effect is due to the inclusion of small scale CMB datasets 
[iniinilTIllinilllllll]. Note from Fig. 4 how our results on WMAP3 + SDSS agree 
with previous analysis with the same data sets [3ll [32] . By comparing our results to 
Fig. 1 of [33 where WMAPl results were reported, we see that models with a natural 
exit from inflation (the region £2 > 0) are now preferred by WMAP3. 

Our analysis tightens the previous bounds on £1 [6l [37]. By combining 
WMAP3-hCMBsmaU with 2dF02 (SDSS) the 2a cl constraint (obtained without 
marginalizing on £2) is £1 < 0.017 (£1 < 0.013). These limits translate on the derived 
parameter Riq < 0.14 {Rio < 0.10) obtained without marginalized on ns- The current 
bound on £1 is almost half of the WMAPl constraint, £1 < 0.031 [J. From this result 
we obtain that the bound on the Hubble scale of inflation driven by a KG scalar field 
derived by WMAP3-f CMBsmall+2dF02: by using the values (not fully marginalized) 
in Fig. ([3]) we obtain H/nip\ < 1.1 x 10^^ which improves the values previously found 
[6l[37] (and implies (V;'"^)^/'* < 2.4 x lO^^ GeV). 

The impact of WMAP3 on simple monomial chaotic models with V{(f)) oc 0" 
strengthens the trend already present in WMAPl. For these models, the HFFs are 
given in terms of the number of e- folds to the end inflation AN as: £1 ~ n£2/4 ~ 

X Our result with WMAPl slightly differs from the value ei < 0.032 6 , since we include also B2K 
in the small scale data set, differently from ^6j. 
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Figure 2. One dimensional marginalized probabilities for cosmological 
parameters obtained by using WMAP3+CMBsmall (black), WMAP3+CMBsmall 
plus 2dF02 (red) or plus SDSS (green). 
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Table 1. Mean values and Icr constraints for the 7 basic parameters and other 
5 derived quantities from the first and three year WMAP release combined with 
the two release of 2dF plus CMBsmall. For ei and Rm the 2ct upper bounds are 
given. 



n/{AAN). The quartic model is now disfavoured at more than 3cr (although the 
agreement of the amplitude of scalar perturbations in such model is not robust with 
respect to the preheating mechanism [34 ). A massive inflaton is in agreement with 
observations, although it moves in the 2a cl with WMAP3 data. 

Power-law inflation 35J obtained by an exponential potential V{4>) = 
Vq exp(— A0/Afpi) (A — y^2/p with a ^ t^) stay on the line £2 = in the (£2 ,ei) 
plane. From the left panel of Fig. 3 and from the relation £1 — — 1) we obtain 
p > 60 at 2ct cl (compared with p > 53 at la cl of Ref. [6]). An interesting issue 
is the appearance of an upper bound (at 2a cl) starts to appear from the left panel 
of Fig. 3, i.e. p < 500 for WMAP3+CMBsmall+2dF02 (or SDSS): as is clear from 
Fig. 4, also this limit may depend on the inclusion of small scale CMB data and for 
its significativity hold the same reasoning made above for the ris = 1 spectrum. Note 
however that the fully marginalized (on non-zero values of £2) value for £1 reported in 
Table I does not have any non-zero lower bound. This upper limit for p, if confirmed 
by future data, can become important in constraining the number of fields which 
support assisted inflation [31], a particle physics realization of power-law inflation. 
The constraints on the exponent of the potential at 2a cl are 0.063 < A < 0.18. 



3. Relaxing the standard consistency condition 

We now discuss the bounds obtained by relaxing the consistency condition (fT2)) and 
by therefore varying 8 parameters (the 7 of Eq. (fT5|) plus the tensor-to-scalar ratio 
To. 01 at the pivot scale fc* = O.OlMpc"^). From a theoretical point of view, this means 
that we are not specifying the Lagrangian p((/), X) for the inflaton and its kinetic 
term X — — Vp(/)V^0/2, but we are restricting ourselves to inflationary spectra at 
first order in HFF. In this way the consistency condition in Eq. (|12p is modified by 
the speed of sound cs for the inflaton [35] : 

s -fi, 2 dp 

— — = -Scsnr , with cg ^ — , (14) 

ros ci^ 
whereas the spectral indices in Eqs. (|6l9p are not modified by cs at first order in 
HFF for constant cs [iQ] . Our analysis with a tensor-to-scalar ratio unrelated to the 
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Figure 4. Constraints by using WMAP3+SDSS (black) vs 

WMAP3+CMBsmall+SDSS (red) on (e2,ei) and (n^ - 1,-Rio) planes at 
1(T and 2(7 level. The right panel shows how the inclusion of small scale CMB 
data drives the spectral scalar index towards red values, disfavouring at just 2a 
cl the Hs = 1 line. The constraints on ng we obtain from WMAP3+SDSS are in 
complete agreemente with |31l 132) . 



tensor spectral index also applies to inflation driven by more than one scalar field, by 
assuming only curvature perturbations after nucleosynthesis (i.e. that isocurvature 
perturbations generated during inflation have completely transferred to the adiabatic 
mode when fluctuations are initialized in the Einstein-Boltzmann code). 

The two dimensional constraints for relevant parameters are presented in Fig. 
([S]). The constraints on the (ei , 62) plane become looser than in the previous section, 
since ei is now unrelated to the tensor-to-scalar ratio. Without marginalizing on the 
scalar spectral index, we obtain r < 0.29 (or Riq < 0.185) at 2cr cl, a bound weaker 
than the one obtained by imposing the standard consistency condition in Eq. (fT2|) 
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Figure 5. Two dimensional contours at la and 2a level obtained by 
WMAP3+CMBsmall+2dF02. The plots reported from left to right are (e2ei), 
(rag — l,i?io), (16ei,r), respectively. In the plane (16ei,r), the red line 
describing the consistency condition r = 16ei is also drawn for comparison. 

(r < 0.27 or Riq < 0.14, see previous section). Note that inflation driven by KG scalar 
field (i.e. cs — 1) is fully consistent with present data. 

4. NEC violating Inflation 

Inflation may occur with a violation of the null energy condition (NEC) |25] . i.e. with 
ei < 0. Since ut = — 2ei, a prediction of NEC violating inflation is a blue spectrum 
for gravitational waves. A similar blue spectrum for gravitational waves may occur 
also in scalar-tensor theories and f{R) gravity theories which are NEC violating in the 
Einstein frame conformally related to the original Jordan one. In order to remain with 
3 inflationary parameters, we shall restrict ourselves to a simple model, just changing 
the sign of the kinetic term of standard inflation, which leads to relation r ~ — 16ei 
in the analysis [41]. The prior used is — 10~^ < ei < —0.1. 

The constraints on the (e2 , ei) for ei < are presented in the left panel of Fig. 
([6]). The contours on (e2 , ci) for NEC violated inflation are completely different from 
standard inflation, although the amount of tensor is slightly loosely constrained with 
respect to standard inflation because of the the blue slope for gravitational waves, as 
can be seen in Fig. ([7]). The best flts for standard inflation and for NEC violating 
inflation have the same value for (= 5673) with WMAP3+CMBsmaU+2dF02. 

From the (e2 , ei) plane, it can be seen that an exponential potential (described by 

the £2 = line) is disfavoured at more than 2a cl. A potential V^cj)) ~ exp (^e2(p^ / 
with 0.02 < £2 ^ 0.08 is consistent with the data. 

5. Second Order HFF Results 

We now allow 632 , 2 7^ in Eq. which, according to inflationary predictions, is 
equivalent to take into account terms which are quadratic in HFFs. Inflationary 
models in which the slow-roll condition is not well satisfied can be described 
analitycally with better accuracy by including the spectra at second order in HFFs 

Let us first try to predict what may happen by inserting the PSs in Eq. Q which 
allow running with respect to the pure power-law PSs analyzed in Section III: at the 
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Figure 6. WMAP3 (red) vs WMAPl (black) la and 2a contours on (e2,ei) 
(left) and {Rio , — 1) for NEC violating inflation. The other data used are 
CMBsmall and 2dF02. 
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Figure 7. Standard inflation (red) vs NEC violating inflation (black) on the 
(ng — 1 , Rio) at la and 2a level. 

price of one additional parameter, the input scalar and tensor PSs are now paraboles 
in ln(fc/fc*), and not straight lines. It is expected that the 7 parameters used in Section 
III change their mean value and broaden their variance in presence of the running, the 
additional parameter. Having in mind the degeneracy in cosmological parameters, it 
is also possible that the best-fit model with runnings may be fairly different from the 
best-fit model without runnings. Within our analysis, we shall see that both these 
possibilities indeed occur with different LSS data sets. 

WMAP Syr analysis included the running of the scalar spectral index as an 
additional parameter [9], neither enforcing a self-consistent tensor-to-scalar ratio to 
second order in HFFs nor including a running for the tensor spectral index. See also 
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Figure 8. One dimensional mean likelihhods (dashed) and marginalized 
probabilities (solid) for ei , £2 i f 3 using a prior —0.2 < es < 0.2 and the GFM. 
We show for comparison results from WMAPl (red) and WMAP3 (black) plus 
CMBsmall+2dF02. 

other studies of the running of the spectral index for WMAPl [42] and WMAP3 [43] . 
We show in Fig. ^ the mean hkchhoods and marginahzed probabihties for 
,^2,^3 with a prior [—0.2,0.2] on €3 obtained by the two WMAP releases. The 
results obtained arc qualitatively in agreement with [B] (for WMAPl) and [37] (for 
WMAP3), although the sets of auxiliary data used are not the same. It is clear from 
Fig. ^ and Eqs. ^ that such prior constrains \as\ < 0.01, although there is a 
tendency to favour positive values for €3, in particular for WMAP3 (as also pointed 
out in [37]). 

In order to have a non trivial result for 63 it is crucial to consider a broader prior 
on this parameter. The justification for this broader prior for £3 is its appearance 
in Eqs. (3-10) not at linear level, but always multiplied by £2- The second order 
formalism is limited by £162 , £2^3 << 1 and not by << 1 (see also f2T for a similar 
claim) . By repeating the same analysis with a much broader prior [§[, we find more 
definite preferred values around £3 ~ 1, as can be seen from Fig. ([9]). In the right 
panel of Fig. ^ we show how the statistical evidence for running has increased from 
WMAPl to WMAP3, although remaining still weak. 

The dependence of ln(P(fc)/Po(fc*)) in Eq. ([3]) and of the coefficients &s,Ti in Eqs. 
(4-10) on de2/dN (— £263) is the reason for a very inefficient exploration of runnings 
of the spectral indexes by the MCMC with basic parameter 63. We find much more 
efficient sampling directly on dt2/dN = 6263 rather than on £3. We present results on 
chains obtained sampling on £263 with prior [—0.5,0.5] in Table 2 and Figs. 9-10. 

At 1(7 we obtain the fully marginalized value of the running of the scalar 
spectral index as = -0.076t|^ S5g (-0-082io ojs) WMAP3-|-CMBsmall-|-2dF05 
(SDSS). The 2cr constraint is -0.141 < ag < -0.013 (-0.142 < as < -0.028) for 
WMAP3+CMBsmaU-t-2dF05 (SDSS). We note that the best-fit cosmological models 
obtained by 2dF05 and SDSS are different, enlarging the difference between the two 
LSS data sets already seen in section III: SDSS prefers a very low value for Hq (as 
also reported in [44]) and the best-fit cosmological model with running is consistent 
with the one without running only at 2a. The cosmological models with and without 
running for 2dF05 are consistent at la: allowing running for the spectral indexes the 
constraint on r < 0.26 for 2dF05 is significantly related to r < 0.50 at 2a cl. 

Since £2 £3 enters only in the scalar running, the magnitude of the running of the 
tensor spectral index is constrained at posteriori to be larj ~ 0(10"'^ — 10^*), as is 
also clear from Eqs. ([9]) and from Fig. (fTO)) . 

Although our combination of data sets and our method of constraining 

§ We have used [-1, 2] for WMAP3+SDSS and [-2, 5] for WMAPl-l-SDSS and WMAP3+2dF. 
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Figure 9. First three panels: one dimensional mean (dashed) and marginalized 
(solid) likelihoods for ei,€2,e3 (left, middle and right panel, respectively) 
obtained by using the two different analytic approximation for the inflationary 
power spectra with WMAP3 + CMBsmall + SDSS: GFM (black) and MCE 
(red). The last panel on the right shows the two dimensional constraints on 
the (ns - 1, as) plane by using the GFM for WMAPl (black) and WMAP3 (red) 
plus CMBsmall and SDSS. 
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Table 2. Mean values and la constraints for the 8 basic parameters and other 5 
derived quantities from WMAP three year data + CMBsmall + 2dF05. Note that 
for ei and Rio the 2a upper bounds are given. The 2a constraints on running of 
the spectral index is; —0.141 < ag < —0.013. 



inflationary parameters differ from those used by the WMAP team [9 , our resuhs 
are in agreement with those reported in [9] for WMAP3 plus other small scale CMB 
data sets. 

6. Conclusions 

We have analyzed the observational constraints from CMB and LSS on the inflationary 
expansion. By parametrizing the inflationary spectra as power-law (i.e. to first order 
in slow-roll parameters), we have updated the WMAPl constraints of Ref. [6]: the 
better determination of reionization obtained by the 3 year data release disfavour now 
both the Harrizon-Zeldovich model (with no tensor) as inflationary models predicting 
ns — I with respect to the first year data. Our analysis shows that models with 
a natural exit from inflation (which stay in the £2 > region of the (e2 ,ei) plane 
according to [33]) are now favoured by WMAP3. 

When running is included, we have shown how the prior on £3 plays a crucial 
role in re-obtaining the WMAP team results within an analysis employing analytic 
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Figure 10. One dimensional marginalized likelihoods for 8 basic parameters and 
7 derived parameters obtained by the GFM formula. We show for comparison the 
constraints by using WMAP3+CMBsmall plus 2dF05 (black) or plus SDSS (red). 




Figure 11. Constraints by using SDSS (red) vs 2dF05 (black) plus 
WMAP3+CMBsmall on (from left to right and top to bottom) (e2,ei), 
(ei , dt2/dN) (ris ~ 1 , Rio), (ris - 1 , as), (as , Rio), (os ,at), planes at Ict and 
2a level. 
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power-spectra depending on (H , e^) with i = 1 — 3. Only allowing a broad prior on €3 
a fairly large negative running in agreement with the WMAP results is obtained. We 
have found that is more convenient to use de2/dN = €263 rather than £3. The value 
for the running as found here implies as ln(A;/fc*)/2 > (ng — 1) for Fourier modes 
far from the pivot scale. It would be therefore interesting to perform a comparison of 
data with inflationary predictions at third order in the HFF. Note that the value of 
the running found here is different from zero at 2a level, and its constraints may 
also be affected by Lya data [HI [46], which we do not consider here. 

We have also studied the impact of theoretical priors on constraining the 
inflationary expansion: still little is known about high-energy physics and prejudices 
in interpreting cosmological data may hide interesting physics. 

We have first relaxed the hypothesis of a canonical scalar field as the inflaton, 
which leads to the consistency relation in Eq. (fT2|) . By considering a generic 
Lagrangian, the constraints on the (ei , €2) plane degrade significantly, while the limits 
on r depend on the slope of the tensor spectrum. This analysis shows that a KG 
inflaton (with sound speed cs = 1) is well inside the confidence contours of the present 
constraints, although the present data do not allow to constrain cs- 

We have then explored how CMB and LSS constrain inflationary models with a 
blue spectrum, a more radical departure than considering r as a free parameter. We 
have restricted ourselves to KG field with the "wrong" sign for the kinetic term, in 
order to do not consider the tensor-to-scalar ratio as an independent parameter. We 
find that both standard inflation and NEC violating inflation are equally good fits to 
the data (the best fits for these two classes of models have the same comparing with 
WMAP3-hCMBsmaU-h2dF02). For nr > 0, the constrained region on the (ei , 62) are 
different from the standard case, leading to 0.02 < £2 ^ 0.08 with |£i| < 0.3 at 2a cl. 
A simple exponential potential potential is disfavoured at 2 cr cl and larger values for 
r are allowed with respect to the standard case with < 0. 
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